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In their recent paper "Faraday rotation of the cosmic microwave background polarization by a 
stochastic magnetic field", Kosowsky et al. |2j have commented about our paper in which we 
derived very strong limits on the amplitude of a primordial magnetic field from gravitational wave 
production. They argue that our limits are erroneous. In this short comment we defend our result. 
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In Ref. Q] we have shown that, if a magnetic field 
is present on super-horizon scales in the early universe, 
its power is very efficiently converted into gravitational 
waves during its evolution from super- to sub-horizon 
scales. We used this fact to derive stringent limits on 
the amplitude of a magnetic field created before the nu- 
cleosynthesis epoch. 

In their recent paper Kosowsky et al. state that 
our limits are not valid. In the discussion section, they 
argue "... the expansion rate of the universe is the same 
whether energy density is converted from magnetic fields 
into gravity waves or not, since the energy density of 
both scale the same way with the expansion of the uni- 
verse. So the actual constraint is on the total radiation 
energy density in the magnetic field, which is constrained 
to be about 1% of the total energy density in the usual 
manner... The corresponding limit on the total comoving 
mean magnetic field strength is around 10 _8 Gauss, not 
the 10~ 27 Gauss claimed in Q." 

We now explain why this conclusion is wrong. We em- 
ploy the same notation convention as 0, and we always 
consider the comoving amplitude of the magnetic field. 
For magnetic fields with spectral index n > —3, the mag- 
netic field energy at wave number k is given by 
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where B\ is the magnetic field amplitude at some fixed 
reference scale A. This energy spectrum is always blue, 
and therefore dominated by its value at the upper cutoff 
k c . This cutoff scale is time dependent, k c (rj). We set the 
magnetic field to zero on scales which are already sub- 
horizon at the time 77* of formation of the magnetic field, 
since we cannot be sure that its spectrum is a power law 
on these very small scales. Therefore, the upper cutoff 
at the time of formation of the magnetic field is given 
by fc* = 777 1 (where rj denotes conformal time). This 
assumption is a conservative one for the derivation of 
our result. At later times, the magnetic field is damped 
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on scales smaller than a time dependent damping scale, 
which gives us the cutoff kr> (77) 0,13 ■ We therefore obtain 
the cutoff function 



fc c (77) = mm(fe,,,£; D (r7)) 



(2) 



Of course at formation fc* = I/77* -C kr)(r}*), while at 
later time feu (77) is decreasing, and eventually becomes 
smaller that /c*. The magnetic field energy density at a 
given time 77 is therefore given by 
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In our paper Q , we have shown that at the time a given 
scale crosses the horizon, and for the maximally allowed 
magnetic fields (Q B ~ ^ rad , see Eqs. (24) and (26) in 
Ref. Q), a considerable fraction 1 of the magnetic field 
energy density is converted into gravitational wave en- 
ergy density: 
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As time goes on, the magnetic field is damped on sub- 
horizon scales k > fee (77) ^> 1/rj, while the gravitational 
waves are not damped, since after formation they no 
longer interact with the matter and the radiation in the 
universe. This is the main point which Kosowsky et al. 
have missed. The magnetic field density parameter at 
nucleosynthesis is given by 
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1 For some values of the spectral index we obtain more energy 
in gravity waves than in the magnetic field. This comes from 
the fact that we linearize the problem and therefore do not take 
into account backreaction. We expect the correct fraction of the 
energy in gravity waves to lie between 50% and 100% of the 
magnetic field energy density. 
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where we have integrated up to the cutoff at nucleosyn- 
thesis, A:D(77 nuc ); while the gravitational wave density pa- 
rameter is 
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where we integrate up to the cutoff corresponding to the 
time of creation of the magnetic field: gravitational wave 
production for a magnetic mode k takes place before hori- 
zon crossing, before the magnetic field is damped by in- 
teraction with the cosmic plasma. A considerable part of 
the magnetic energy is therefore "stored" in gravitational 
waves. 

From (0J and 0, one can see that the ratio be- 
tween the two energy densities is ris(?7nuc)/^G — 
(&D(?7nuc)A*) n+3 - Only for n ~ -3 or 77* ~ r] nuc this 
factor is of order unity; in this case, one can apply the 
nucleosynthesis bound indifferently to SIq or £Ib , and one 
gets the same constraint. In most cases instead, the ratio 
is rather huge. Let us consider the example of magnetic 
field generation at the electroweak phase transition. In pj 
we calculate k£>(r] DUC ) ~ 6 x 10~ 7 sec _1 ~ 60 pc _1 , and 
rj cw ~ 4 x 10 4 sec. Taking into account that electroweak 
magnetic field generation is causal (not inflationary), and 
therefore n=2 (see 0), we obtain 
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and by no means one! If the magnetic field is generated 
during inflation, one is no longer forced to have n = 2, 
but can have arbitrary values of n > — 3. If we take 
n ~ 0, for an inflation scale of 10 15 GeV, we have fe* = 
l/ryint ~ 10 13 /r7 ow , and we obtain 
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We can conclude that, applying the nucleosynthesis 
bound on fig, we find much stronger constraints on B\, 
due to the fact that a seizable fraction of the magnetic 
field energy is converted into gravitational waves before 
the damping process. 

Apart from not taking into account this damping, there 
is a second point that has been missed in Kosowsky et 
al. We specifically pronounce a limit for the amplitude 
of the stochastic magnetic field smoothed over a scale 
A ~ 0.1 Mpc. On the contrary, they talk about 'the total 
comoving mean magnetic field', which is largely domi- 
nated by its value on small scales, hence B(ko(j]nuc))- 



The value of the field at this scale is limited to a few 
10 _8 Gauss by the constraint Qb < 0.1fi ra d at nucle- 
osynthesis. But this field value has no relevance, for two 
reasons. First of all, because the damping scale will 
grow, which means that B(ko(j]nuc)) will be damped 
away before it can ever give rise to magnetic fields in 
galaxies. Typically, the highest mode which survives 
damping is kn{rj rec ) ~ 10 Mpc -1 , much smaller than 
&D(?7nuc) H, 0| • Secondly, the scale relevant for mag- 
netic fields in galaxies and clusters is A ~ 0.1 — 1 Mpc, 
and we have thus formulated limits for this scale. If 
5(Mw))Mw) 3/2 ee B kD < 10" 8 Gauss, the limit 
on the scale A 3> l/fco^nuc) is much smaller, namely 
B X = B(k = 1/A)A- 3 / 2 = 5 feD (Mw)A)- (n+3)/2 < 
10~ 8 Gauss x io- 6 ("+ 3 )/ 2 . For a spectral index n = 2, 
for example, B\ is smaller than the maximal field by a 
factor of ~ 10 15 , namely B\ < 10~ 23 Gauss. 

In conclusion, even if the magnetic field and gravita- 
tional wave energy densities scale in the same way with 
the expansion if the universe, applying the nucleosyn- 
thesis bound on the induced gravitational wave energy 
density gives a much stronger constraint on the ampli- 
tude of the magnetic field. Both the magnetic field and 
gravitational wave energy spectra are blue, and therefore 
dominated by their value at the upper cutoff. However, 
the upper cutoff for the gravitational wave spectrum is 
much higher than the one for the magnetic field spectrum 
at the epoch of nucleosynthesis: r/^ 1 ^> fczj^nuc)- The 
reason for this being, that the conversion of magnetic 
field energy into gravitational wave energy takes place 
when a given mode enters the horizon, before the mag- 
netic field is dissipated by interaction with the cosmic 
fluid. 

Furthermore, the interesting limit is not the one on the 
'mean magnetic field' which is dominated by the value at 
the smallest scale, but the limit on the field amplitude 
at some scale A which is relevant for galactic magnetic 
fields, and certainly has to be larger than the damping 
scale at the redshift of galaxy formation. 

Of course, these limits apply for magnetic fields gener- 
ated before the epoch of nucleosynthesis. 
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